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EDITOR'S  SUMMARY 

In  this  estimation  of  impact  pressures  of  a  drop  with  a  solid  surface,  the 
deformation  of  the  surface  of  the  drop  and  the  onset  of  post-collisional  motion 
of  the  liquid  are  considered.  Equations  of  maximum  pressure  are  developed  and 
their  dependence  on  the  extent  of  wetting  of  the  solid  surface  shown  to  be 
significant.  They  are  not,  however,  greatly  affected  by  the  shape  of  the  free 
surface  formed  by  the  drop  on  impact.  It  is  also  shown  how  excessive 
flattening  of  the  drop  leads  to  a  reduction  of  impact  pressures. 
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The  results  set  out  in  this  article  are  a  further  continuation  of  work* 
which  gave  a  simplified  scheme  for  estimating  the  impact  of  a  drop  against  a 
solid  surface. 

In  that  work  the  calculation  was  based  on  the  assumption  that  at  the 
moment  of  impact  the  drop  assumed  the  form  of  a  hemisphere.  In  the  present 
article  that  constraint  is  abandoned  and  the  problem  of  the  impact  of  a  drop 
with  a  surface  of  arbitrary  shape  after  impact  is  considered. 

In  this  case  the  axi symmetrical  problem  is  replaced  by  a  plane  one.  This 
schematisation  (as  shown  in  Ref. 2)  does  not  have  any  significant  influence  on  the 
magnitude  of  the  impulsive  pressures  obtained. 

(1)  At  the  moment  of  impact  of  a  drop  against  a  solid  surface  a  change  in  the 
shape  of  the  surface  of  the  drop  must  occur  as  well  as  motion  of  the  liquid 
within  the  drop.  Instead  of  considering  this  complex  unsteady  process  we  shall 
consider,  as  in  the  study  of  the  impact  of  solid  bodies  ,  the  limiting  motion 
which  occurs  in  step-wise  change  of  velocity  of  a  drop  (in  a  system  of 
coordinates  coupled  with  the  liquid,  its  velocity  before  impact  will  be  every¬ 
where  equal  to  zero) .  So  as  to  take  into  account  the  change  in  the  surface  of 
the  drop  on  impact  we  shall  consider  the  flow  arising  in  the  drop  deformed  by 
impact.  Thus,  the  schematisation  of  the  process  of  impact  of  the  drop  consists 
in  the  fact  that  the  continuous  process  of  onset  of  internal  motion  in  the 
liquid  and  deformation  of  the  drop  is  split  into  two  separate  phases: 

(I)  deformation  of  the  surface  of  the  drop;  (2)  onset  of  post-collisional 
motion  of  the  liquid. 

Assuming  that  the  liquid  is  ideal  and  incompressible  we  find  that  the 

problem  of  calculation  of  the  motion  of  the  liquid  in  the  drop  following  impact 

comes  down  to  determining  the  velocity  potential  V  for  the  following  boundary 
.  .  3 

conditions  :  at  the  free  surface  ^  ■  0,  at  the  solid  wall  the  normal  component 
of  velocity  is  equal  in  magnitude  and  opposite  in  direction  to  the  velocity  of 
the  drop  prior  to  impact. 

In  what  follows  we  shall  consider  plane  flow,  substituting  for  a  drop  an 
infinitely  long  liquid  circular  cylinder  moving  prior  to  impact  in  the  negative 
direction  of  the  axis  Ox  (the  axis  Oy  representing  a  section  of  a  solid 
plane  wall),  see  Fig.l. 
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The  problem  of  calculating  this  flow  will  then  consist  in  determining  the 
complex  flow  potential  w  ■  <p  +  ii|»  ■  potential,  \p  ■  stream  function)  in  the 
region  bounded  by  the  plane  solid  wall  AOB  and  the  contour  of  the  free  surface 
ACB  for  the  following  boundary  conditions: 

(1)  at  AOB  the  velocity  component  along  the  x  axis  v^  «  c 
(c  -  pre-collisional  velocity  of  the  drop); 

(2)  at  ACB  the  potential  ^  ■  0. 

(2)  In  order  to  determine  the  dependence  of  w  on  z  we  shall  map  the 
region  AOBC  on  the  interior  of  a  unit  semicircle  in  the  region  of  an 
additional  variable  £  *  £  +  in  (Fig. 2).  Let  us  first  find  the  additional 
function 


a 


dw  .  dw 
d?  "  15  dz  d? 


At  the  solid  wall  AOB 


dz  .  dz 

dC  "  "  1  d? 


-  iy(£)  • 


(2-1) 


(2-2; 


The  function  u(£)  which  is  positive  everywhere  within  the  segment  AOB 
is  determined  by  the  shape  of  the  contour,  and 


dw 

dz 


(2-3) 


Here  v^  and  are  the  fl°w  velocity  components  along  the  x  and  y  axes. 

Thus 


Re  fl  ■  c£ii(€)  at  AOB  . 


Since  ^  ■  0  on  the  arc  ACB,  on  this  arc 


dw 

dC 


+ 


dw 

d? 


e-i  >rg  {  _ 


n 


±  i 


(2-4) 


and  hence 


Re  0 


0  . 


(2-5) 


-T 
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Making  use  of  this  we  can  continue  the  function  ft  analytically  through 
the  arc  ACB  on  to  the  whole  upper  half-plane  of  the  region  of  change  of  C. 
After  this  we  get  for  Re  ft  in  the  interval  (-1,  «,  +1)  of  the  axis  0£ 


Re  ft(0  -  -  Re  ft(£_1) 


f  vU~') 


(2-6) 


Knowing  the  value  of  Re  ft  over  the  whole  of  the  axis  Q£  we  can 
determine  ft  by  means  of  the  Schwartz  formula^ 


-1  ,  ®  1 
f  i  dT  +  f  w(t~‘)  dT  _  [  aii) 

I  T(T  -  C)  J  T(T  -  C)  ^  J  T  -  C 

1  -1 


dt>  +  ic  (2-7) 


(C  is  n  real  constant  which  requires  to  be  determined).  Substituting  for  ft 
the  expression  (2-1)  we  get 


dw 

dz 


_e_ 

it 


■ 


E 


T " 5  5(t  -  rS 


y(i)dT 


>*?  • 


(2-8) 


formula 


For  t  ■  5C(-I,  +1)  we  obtain  from  (2-8)  by  means  of  the  Sokhotski 
4 


dp 

d? 


c 

TTC 


1  r- 


/ 

-1 


tu(t)  - 


T  -  5 


M(T)  1 

•'  -  'J 


dT  +  KVU)  In 


*  ~  •  (2-9) 


Since 


d£ 

d? 


d£ 

dy 


dz 

dt 


on  AOB 


& 

dy 


0  at  the  point  0 


then 


d<fi 

d5 


0  when  £  ■  0  . 


From  this  it  follows  that  C  ■  0. 

Assuming  that  w  -  0  at  point  A  and  integrating  (2-8)  we  obtain  for  the 
complex  potential  at  an  arbitrary  point  within  the  liquid 
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p(T)dTd£  = 

If  the  map  of  the  region  (z)  on  the  semicircle  (?)  is  known  then  the 
expression  (2-10)  together  with  the  function  z  =  z(?)  will  give  the  dependence 
of  w  on  z  in  parametric  form. 

(3)  The  impulsive  pressure  on  a  solid  surface  can  be  determined  from  the 
formula'* 


-5/* 

-1 


(t)  In 


i  -  ;t 

t  -  ; 


dr.  (2-10) 


p  =  - 


(3-1) 


Passing  to  the  limit  in  (2-10)  when  ?  ->•  £  we  obtain  for  the  complex 
potential  at  a  solid  surface 


w 


m(t)  In  *  ^  -  dt  +  ic(y  -  a) 


(3-2) 


Substituting  in  (3-1)  the  value  of  the  potential  'f  taken  from  this 
formula  and  the  value  of  u(x)  from  (2-2)  we  define  the  impulsive  pressure  at 
any  point  of  a  solid  surface  as 


P 


C  -  T 

1  - 


di  . 


(3-3) 


It  follows  from  this  expression  that  at  the  boundaries  of  a  wetted  surface 
the  impulsive  pressure  is  equal  to  zero.  The  maximum  pressure  which  arises  at 
the  centre  of  the  drop  will  on  the  basis  of  (3-3)  be  equal  (substituting  5  *  0 
in  (3-3)  and  carrying  out  partial  integration)  to  / 


P* 


2pc 

7T 


H  d> 


(3-4) 


(4)  Let  us  consider  how  the  magnitude  of  the  maximum  impulsive  pressure 
depends  on  the  shape  of  the  drop  after  impact  and  the  magnitude  of  the  wetted 
surface.  For  this  we  calculate  the  maximum  impulsive  pressure  p^  for  several 
different  contours  of  the  free  surface  following  impact:  (1)  arc  of  a  circle; 
(2)  ellipse;  (3)  inversion  of  ellipse  relative  to  a  circle  passing  through  its 
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minor  semi-axis.  Each  of  these  contours  is  determined  by  two  parameters,  one  of 
which  is  the  half-width  of  the  wetted  surface  a. 

The  second  parameter  is  defined  such  that  the  area  of  section  of  the 
liquid  cylinder  remains  constant  before  and  after  impact  (condition  of 
incompressibility).  As  the  second  parameter  we  shall  take  the  maximum  width  of 
the  drop  in  the  direction  of  the  abscissa  axis  b  (Fig.l).  For  the  magnitude 
of  the  pressure  we  get  the  following  formulae  respectively  (if  the  radius 
of  section  of  the  liquid  cylinder  before  impact  is  denoted  by  R): 

(1)  Arc  of  circle 

2pc  f  1  -  [(1  -  x) / ( 1  +  t)]  dt  4pc 

'  i  !♦[<!-  0/0  *  T)]29"  T  '  * 

....  (4-1) 


/  1  -  t29/* 

I  (1  +  t26/lT)(l  -  t2) 


Here  6  is  the  angle  between  a  tangent  to  the  arc  and  the  solid  surface  at 
their  point  of  intersection  measured  from  within  the  liquid  (boundary  angle). 
Its  relation  to  a  is  given  by 

a  ■  K^aineJw^£r I#  •  <4-2) 

(2)  Ellipse  (a  and  b  ■  major  and  minor  semi-axes  of  the  ellipse 
respectively) 


P* 


fj  C< 


a  +  b)e  9  +  (a 


b)e  ] 


B1  dT  2pc  /a 


cos  B 


dt  (4-3) 


K 


(4-4) 


.  /(I  ♦  x: 


)2k2  -  1  \ 
T2)  ) 


B  - 


y  F  ^k,  ar 

7TT  F  fk' ,  arc  sin  — -  \ 

,W  \  k'(l  +  t  nj 


when  t  <  x  (4-5) 


when  t  >  x  .  (4-6) 
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Here  K  and  F  are  respectively  the  complete  and  incomplete  elliptic  integrals 
of  the  first  kind,  and  k  and  k'  the  module  and  additional  modulus  of  this 
integral,  the  modulus  k  being  determined  from  the  major  semi-axis  of  the 
ellipse  by  the  relationship 


K(k') 

K(k) 


(4-7) 


We  should  observe  that  even  when  (a//2  R)  >  1.5  the  difference  between 
(1  -  k')/k  and  1  is  very  small  and  k'  <3  !.  Thus  in  (4-3)  the  second 
integral  can  be  discarded  and  the  expression  for  B  simplified. 


Then 


2pc 


/  2  .2 

/a  +  b 


l-k' 


cos 


2K(k) 


it  „  /,  2ik'  \ 

F  Ik,  arc  cos  - 


(4-8) 


In  this  case 
k'  -  4e-R(k> 


K  -  7T 


2  In 


2  2' 

aT  +  2R 

2  2 

a  -  2R 


hi 


K(k')  -  1  when  -i- >  2 
2  /2R 


(3)  Inversion  of  ellipse 


dt 


(b  -  a)t  +  (b  +  a) 


4cpab 


TT  A2  - 


arctg 


b  + 


(4-9) 


where 


b  = 


(4-10) 


When  a  -  b  all  the  contours  considered  become  a  semicircle  of  radius  /2  R 
and  for  the  maximum  pressure  we  get  the  value  given  in* 


P** 


2pca  _  2  /l  Rpc 

TT  IT 


(4-1!) 


(5)  The  dependence  of  the  maximum  pressure  on  the  width  of  the  wetted  part  of 
the  solid  surface  for  the  types  of  contours  considered  is  shown  in  Fig. 3  where 
the  graphs  for  the  change  of  the  relation  o  *  P*/p**  are  shown  as  a  function 
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of  a  ■  a//2  R.  These  graphs  show  that  the  magnitude  of  the  impulsive  pressure 
is  considerably  influenced  by  the  magnitude  of  the  wetted  solid  surface;  the 
shape  of  the  free  surface  formed  by  the  drop  on  impact  does  not  greatly  affect 
the  magnitude  of  the  impulsive  pressure.  At  the  same  time  the  results  of  these 
calculations  show  that  excessive  flattening  of  the  drop  and  also  insufficient 
flattening  of  it  lead  to  reduction  of  the  impact  pressures.  The  magnitude  of 
as  determined  from  (4— 11)  gives  the  value  of  the  upper  limit  of  the 
impulsive  pressures  arising  on  impact  of  the  drop. 
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